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Abstract 

We describe exchange-only universal quantum computation and leak- 
age reduction in the 3-qubit decoherence free subsystem (DFS). We dis- 
cuss the angular momentum structure of the DFS, the proper forms for 
the DFS CNOT and leakage reduction operators in the total angular mo- 
mentum basis, and new exchange-only pulse sequences for the CNOT and 
leakage reduction operators. Our new DFS CNOT sequence requires 22 
pulses in 13 time steps. The DFS leakage reduction sequence, the first 
explicit leakage reduction sequence of its kind, requires 30 pulses in 20 
time steps. Although the search for sequences was performed numerically 
using a genetic algorithm, the solutions presented here are exact, with 
closed-form expressions. 

1 Introduction 

Interest and research in semiconductor quantum dots for quantum information 
processing has continued to grow since the original proposal by Loss and Di- 
Vincenzo [IJ[2]. In semiconductor quantum dot systems, with a single electron 
spin as a qubit, two qubit gates using electron exchange interactions can be per- 
formed on sub- nanosecond time scales [3]. In contrast, single qubit operations 
based on electron spin resonance may be two orders of magnitude slower or 
more [H [5] . The slow and technically challenging single electron operations can 
be avoided by using an encoding in which exchange interactions give encoded 
universality. The possibility of universal computation using only the exchange 
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interaction, and its connection to decoherence free (DF) subsystems, has been 
demonstrated in [6j[7j|8]. Bacon et al. [6] describe the 4-qubit DF subspace, 
giving exchange-based Hamiltonians that generate encoded universal computa- 
tion. Kempe et al. |S] give a general theory for decoherence free subspaces 
and subsystems, and prove the universality properties of these encodings em- 
ploying only exchange or other two-body interactions. Experiments on a single 
three-spin exchange-only encoded qubit have also recently been reported [9]. 

Any implementation of exchange-only quantum computing, however, re- 
quires explicit gate sequences for a universal set of encoded gates. For the 
3-qubit decoherence free subspace, DiVincenzo et al. [TU] found explicit ex- 
change gate pulse sequences for single encoded qubit gates, as well as for an 
encoded two-qubit gate locally equivalent to a CNOT. (A two-qubit gate is lo- 
cally equivalent to a second two-qubit gate if they differ only by single qubit 
gates [TT].) Kawano et al. [H] showed that DiVincenzo's numerically obtained 
locally equivalent CNOT sequence using 19 gate pulses in 13 time steps in fact 
approximates a true, analytical, solution. 

The DiVincenzo 19 gate locally equivalent CNOT pulse sequence is valid 
only for the 3-qubit DF subspace, and not for the entire 3-qubit DF subsystem. 
Computation in the DF subsystem offers two major advantages compared to the 
subspace: initialization of encoded subsystem states requires no magnetic field, 
and the subsystem states are immune to all global decoherence mechanisms, not 
just global decoherence in a single direction. To exploit these features of the 3- 
qubit subsystem, an analytic CNOT gate sequence for the entire subsystem was 
found by Bonesteel et al. [13], which requires approximately 50 exchange gates. 
Using a genetic algorithm, we have significantly improved on the Bonesteel et 
al. solution, finding a new analytic pulse sequence for the encoded subsystem 
CNOT that requires just 22 exchange gates in 13 time steps. Additionally, our 
solution is for the full encoded CNOT, and not just the locally equivalent CNOT. 
The 4-qubit DF subspace possesses the same magnetic field-free initialization 
and global decoherence immunity as the 3-qubit subsystem. An encoded CNOT 
gate sequence has also been found for the 4-qubit encoding |14j . requiring 50 
exchange gates in 27 time steps. With a CNOT gate sequence that is more than 
twice as expensive as our 3-qubit CNOT solution, and the additional overhead 
of an extra physical qubit per encoded qubit, the 4-qubit encoding is apparently 
inferior to the 3-qubit encoding. 

Computation with the 3-qubit DF subsystem requires not only a univer- 
sal set of encoded gates, but also an effective means of recovering from errors. 
Though the subsystem protects against global decoherence, local decoherence 
mechanisms such as nuclear hyperfine and electron-electron dipole-dipole cou- 
pling still exist, and give rise to both encoded errors and leakage from the 
encoded subsystem. Encoded errors can be corrected using standard quantum 
error correction procedures [15], but leakage errors must be converted, or re- 
duced, to encoded errors in the course of applying quantum error correction. 
The incorporation of leakage reduction units into fault tolerant quantum error 
correction circuits is described in [16]. Kempe et al. [8] describe a "SWAP If 
Leaked" (SIL) operator that performs leakage reduction for the 3-qubit DF sub- 
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Table 1: Quantum numbers of commuting operators that uniquely specify all 
basis vectors in a three qubit, eight-dimensional Hilbert space. S is the total 
spin of all three qubits and specifies whether the DFS qubit has leaked; S = i 
is unleaked, S = ^ is leaked. Si 2 is the total spin of the first two qubits and 
gives the logical or encoded state for unleaked states. S z is the total spin-z of 
all three qubits and is the gauge quantum number. The top line gives the index 
labels of the basis vectors, which correspond to Eqs. {I])-©. 



space with only exchange gates: given a potentially leaked encoded qubit A and 
a fiducial, unleaked encoded qubit B, the SIL operator leaves A unchanged if 
it is unleaked, and replaces A with a valid encoded state if it is leaked. Kempe 
et al. discuss the action of an SIL operator on encoded basis states but do not 
give an explicit pulse sequence for SIL. Here we present the first exchange-only 
pulse sequence for leakage reduction in the 3-qubit DF subsystem, requiring 
30 exchange gates in 20 time steps, with closed-form pulse timings. Together, 
our leakage reduction and CNOT pulse sequences provide two of the required 
operations for computation using the 3-qubit DFS. 

The paper is organized as follows. In Section 2 we review the angular mo- 
mentum structure and basis states of one and two 3-qubit DF subsystems and 
define the notation to be used in the following sections. In Section 3 we give the 
form for the encoded CNOT operator in the total angular momentum basis and 
the exchange gate sequence satisfying this form. Section 4 gives the analogous 
form and sequence for the leakage reduction operator. Section 5 describes the 
genetic algorithm used in the numerical searches for the pulse sequences. We 
conclude in Section 6. In the remainder of the paper we use the acronym "DFS" 
to denote decoherence free subsystem. 

2 Angular Momentum Structure of DFS 

The basis states of the Hilbert space containing three qubits can be described by 
three angular momentum quantum numbers S, Si, 2, and S z [Z|- S is the total 
spin of the three qubits, corresponding to the eigenvalues of the operator S 2 = 
S 2 + Sy + S 2 , where Sj = l/2(ay + uj 2 ^ + uj 3 ^) is the total angular momentum 
in the j direction, is the j th Pauli matrix on the n th qubit, and h = 1. The 
eigenvalues of S 2 are S(S + 1). Si ,2 is similarly the total spin of the first two 
qubits, and S z the total z spin of all three qubits. The eigenstates described by 
these quantum numbers are shown in Table [TJ Their corresponding eigenvectors 
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can be written in the computational basis via Clebsch-Gordan coefficients: 



!>= ^(|010) - |100)) =\S )\0) (1) 

2) = ^(|0H) -|101)) = |Sb> |1> (2) 

3) = ^/f |001) |010) |100) =-^(V2|T + )|l)-|T )|0)) (3) 

4) = fa |011) + ^ |101) -yj* |110) =-^(|T )|l)-v^|T_)|0) (4) 

5) = |000) =|I+)|0) (5) 

6) = ^(|001) + |010) + |100)) =i=(|T + )|l)+V2|T )|0)) (6) 

7) = ^(|011) + |101) + |110)) =i(V2|T )|l) + |T_)|0)) (7) 

8) = |111) =|T_)|1) (8) 
where the singlet |So) and triplet |T M ) states are defined as: 

\So) = ^=(|01) - |10>) (9) 

\T+) = |00) (10) 

\To) = ^(|01) + |10)) (11) 

|T_) = |11) (12) 

The singlet and triplets are eigenstates of the 2 operator. 

The 3-qubit DFS is spanned by the first four eigenstates in Table [TJ The 
total spin quantum number S distinguishes between states in the valid unleaked 
subspace (S = |) and leaked states (S = |). Within the unleaked subspace 
a valid state of the 3-qubit DF subsystem has encoded (or logical) quantum 
number Si t 2 and gauge quantum number S z that are unentangled: valid states 
are factorizable states of the abstract (rather than physical) subsystems corre- 
sponding to Si 2 an d S z quantum numbers. Valid DFS subsystem states are 
given by a(j \1) + S |2)) + /3(7 |3) + S |4)). In terms of the Si,2 and S z quantum 
numbers we see the factorizability explicitly: 

a( 7 |l)+*|2))+0( 7 |3)+*|4)) 

= a(rt \S = 1/2, Si j2 = 0, S z = 1/2) + 5 |l/2, 0, -1/2)) 

+/3( 7 11/2,1,1/2) +6 |l/2,l,-l/2)) (13) 

Global decoherence mechanisms couple only to the S z quantum number, and 
thus modify the gauge state only and not the encoded information. A 3-qubit 
DFS state can be initialized as a singlet in the first two qubits and an arbitrary 
spin in the third qubit, giving the state 7 |1) + 6 \2). The initialized state has 
encoded quantum number Si,2 = and an undefined gauge state unentangled 
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with the encoded quantum number. Assuming that the 3-qubit DFS is not 
leaked, measurement of the total spin £1,2 on the first two qubits distinguishes 
between singlet and triplet states and gives the logical content of the DFS. 

Encoded operations on the 3-qubit DFS are performed using the exchange 
interaction between its constituent qubits. The exchange interaction between 

qubits m and n is defined as H^ n = 1/4 (tr^ ■ + ay ■ a y + a z • a z )■ 
The exchange interaction generates the SWAP operation between qubits m and 
n, with a partial swap operation given by 

U%n(p) = e M-™pHm,n); ( 14 ) 

Um,n(P = il) gives the usual SWAP operation up to a global phase. Because 
the exchange interaction generates swaps between qubits making up the 3-qubit 
DFS, it cannot change any total angular momentum quantum numbers: it com- 
mutes with both S and S z (as well as S x and S y ). The exchange interaction 
does change the S\^ quantum number and hence the encoded state: exchange 
between qubits 1 and 2 generates an encoded z rotation, and exchange between 
qubits 2 and 3 generates an encoded rotation about the h = {\/3/2, 0, — 1/2} 
axis. 

The fact that the exchange interaction cannot change total angular momen- 
tum quantum numbers motivates the choice for the set of quantum numbers 
describing two 3-qubit DFS's or six physical qubits. For DFS qubits A and B, a 
valid basis set is simply the "product basis" of the A and B eigenvectors given in 
Tabic [T] None of the quantum numbers of the product basis commutes with all 
of the possible exchange operations between the six physical qubits. However, if 
we use the total angular momentum basis consisting of quantum numbers Stot) 
S z , tot, Sa, Sb, Sa,i.2, and Sb.i,2, the conservation of Stot and S z> tot leads to 
a partial diagonalization of the exchange-constructed operator and the block 
structure described below. Stot is the total spin of all six physical qubits and 
S z , tot is the total z spin of all six physical qubits. Sa {Sb) is the total spin of 
DFS qubit A (B), with Sa = § and Sb = \ corresponding to unlcaked states. 
Sa,\,2 (Sb, 1.2) is the spin of the first two qubits of DFS qubit A (B) and gives 
the logical information encoded in DFS qubit A (B). 

The structure of a six-qubit system is given in [5] . It consists of five spin-0 
subspaces, nine spin-1 subspaces, five spin-2 subspaces, and one spin-3 subspace 
[17] . Each total spin subspace is further divided into 5 z ,tot subspaces. Tables [2]- 
HJshow the quantum numbers for the basis states of Stot = 0, 1, 2, respectively 
(Stot = 3 is not needed in the following). Table[2]gives the five-dimensional spin- 
subspace of the six physical qubit system. Table |3] gives the nine-dimensional 
Stot = lj £ Zl tot = ~ 1 subspace; Stot = 1, S Z) tot = or 1 have analogous nine- 
dimensional subspaces. Table 2] gives the five-dimensional S to t = 2, S Zl t D t = ~ 2 
subspace; similarly, S to t = 2 and S Zi tot — — 1) 0, 1, or 2 have analogous five- 
dimensional subspaces. The 64-dimensional basis vectors can again be written in 
terms of the computational basis of six spin-1/2 particles using Clebsch-Gordan 
coefficients. 

Because the exchange interaction commutes with Stot and S Zi tot, any op- 
erator constructed from exchange gates cannot mix subspaces with different 



5 



1 2 3 4 5 









Si 

S b 



11113 
2 2 2 2 2 
11113 
2 2 2 2 2 




111 
10 11 



Table 2: Quantum numbers for two DFS qubits in the total spin-0 subspace. 
Stot is the total spin of the six physical qubits, S z ,tot the total spin-z, Sa (Sb) 
the total spin of DFS qubit A (B), Sa, 1.2 (Sb, 1,2) the spin of the first two qubits 
of DFS qubit A (B). The top row is the index label of the basis vectors in the 
total angular momentum basis. Basis vectors 1-4 are valid encoded states; basis 
vector 5 is a leaked state. 

Stot an d S Zl tot- Additionally, because the exchange interaction also commutes 
with S Xi tot and S^tot, exchange-constructed operators must have exactly the 
same matrix entries on subspaces with the same Stot but different S Zj tot- That 
is, the matrix entries of an exchange-constructed operator on the Stot = 1) 
S Z) tot = — 1 subspace are exactly the same on the Stot = 1, <S«,tot = and 
Stot = 1, Sz,tot = 1 subspaces, and similarly for S to t = 2, S z , to t = -2, -1, 0, 1, 
2. When written in the total angular momentum basis, an exchange-constructed 
operator is a block diagonal matrix, whose diagonal blocks are one 5x5 spin-0 
block, three identical 9x9 spin-1 blocks, five identical 5x5 spin-2 blocks, and 
a spin-3 block consisting of a 7x7 identity matrix times a phase factor; none 
of these blocks couples to another. If we define the generator of swap to be 
H™ n = H% >n - 1/4, so that exp(-/7rtf^„) gives a full SWAP, including the 
global phase, the swap-generated group is SU(5) x S/7(9) x St/(5) x £7(1). The 
Stot = 3 subspace transforms as the identity under swaps, so the U(l) phase 
factor must be exp(— i#[St t(St t + 1) — 12]). The algebraic origin of the block 
structure is described in [71 H] ■ Constructing DFS operators on the two DFS 
qubits amounts to finding exchange gate sequences that satisfy the desired forms 
of the block matrix in the total angular momentum basis. 



The CNOT operation on two unleaked 3-qubit DFS's must perform a CNOT 
gate on the logical DFS information contained in the quantum numbers Sa,i,2 
and Sb,i,2, independent of the gauge state. We seek an exchange pulse se- 
quence, i.e., a sequence of unitary operators generated by the exchange inter- 
action, whose product in the total angular momentum basis gives the required 
CNOT operation on the encoded quantum numbers. Because of the structure 
of exchange-constructed matrices described in Section [2] we need only constrain 



3 DFS CNOT Pulse Sequence 
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Table 3: Quantum numbers for two DFS qubits in the total spin-1, S z .tot = — 1 
subspace. Basis vectors 6-9 are valid encoded states; basis vectors are leaked 
states, with one or both of the constituent DFS qubits leaked. Basis vectors 10 
and 11 are unleaked in DFS qubit A. See Table [2] for operator definitions. 
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Table 4: Quantum numbers for two DFS qubits in the total spin-2, S z ,tot = — 2 
subspace. Basis vectors 15 and 16 are unleaked in DFS qubit A. See Table [2] 
for operator definitions. 
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the S tot = 0, S'a.tot = 0, 5 x 5 block and S tot = 1, Ss.tat = -1, 9 X 9 block 
appropriately for the DFS CNOT. In the S to t = 0, S z>tot = block the CNOT 
matrix must take the following form: 
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where 9q and 0o are arbitrary phases, and the rows and columns correspond to 
the basis vectors 1-5 in Table [2j The upper left 4x4 block is the usual CNOT 
operation on the unleaked, logical states; unitarity requires that the leaked state 
with Stot = 0' basis vector 5, be uncoupled from the unleaked states. In the 
three 9x9 spin-1 blocks the CNOT matrix must be 
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Cl,l Cl,2 Cl,3 Cl,4 Cl,5 

C2,l C2,2 C2 : 3 C2,4 C2,5 

C3,l C 3;2 C 3:3 C 3i4 C 3j5 

C4,l C4 ; 2 C4 : 3 C4 : 4 €4,5 

C5,l C5,2 C5,3 C5,4 C5,5 / 



(16) 



where again 0\ is an arbitrary phase and the rows and columns correspond to 
basis vectors 6-14 in Table [3J The upper left 4x4 block must be the CNOT on 
the encoded basis vectors, and the lower right 5x5 block {cjj} is an arbitrary 
unitary matrix on the leaked states 10-14. Again, the leaked and unleaked 
states must not couple. The three 9x9 spin-1 blocks are automatically iden- 
tical when constructed out of exchange gates, so they need not be constrained 
separately. The CNOT matrix on the spin-2 and spin-3 subspaces is completely 
unconstrained, aside from unitarity, which is automatically satisfied with the 
exchange gates. 

An objective function for the CNOT search is constructed from the con- 
straints shown in Eqs. (fl~5|) and (fl~6|) . Let U be the product of exchange gate 
unitaries in the total angular momentum basis, ordered so that basis vectors 
1-5 correspond to Table [5] and basis vectors 6-14 correspond to Table [31 The 
objective function for the CNOT search is 



/cnot(C^) = 



U. 



2,2 



U 3 .4 + U 4 , 3 \ - -| 176, 6 



7,7 ■ 



8,9 



09,8 1 ■ 

(17) 

Since U is unitary by construction, its entries have modulus at most 1. The 
objective function is zero only when Ui t i, c/2,2, ^3,4, and U^^ have modulus 1 
and a common phase, and E/6,67 U-j^ 1 £/s,9> and C/9,8 also have modulus 1 and 
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Figure 1: Twenty-two pulse, 13 time step, exchange gate sequence for DFS 
CNOT. DFS qubits A (control) and B (target) are arranged as shown. Sub- 
scripts on A and B label the constituent physical qubits. Each gate corresponds 
to an exchange unitary U™ n (p) [Eq. ([13)) ] . with the swap powers p displayed 
explicitly in the gates, p — 1 corresponds to a full SWAP operation, up to a 
global phase. p\ = arccos(— l/-\/3)/7r and pi = arcsin(l/3)/7T. 



a (generally different) common phase. This objective function is used in the 
genetic algorithm described in Section [5l The two DFS qubits A (control) and 
B (target) are laid out in a linear array, in the order A3, A2, Ai, B\, B2, -B3. 
Only nearest neighbor exchange gates are permitted. 

Figure Q] shows the best solution found by the genetic algorithm requiring 22 
pulses in 13 time steps. Though the genetic algorithm uses approximate (finite 
precision) numbers, the final solution found is analytic, with the exchange gate 
powers shown in the figure. Writing the exchange unitaries given in Figure [1] in 
the total angular momentum basis and taking their product yields the following 
matrix for the spin-0 subspace: 
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(19) 



The spin-0 and spin-1 subspaces have the same global phase. The matrices in 
Eqs. ([T8l) and ([T9"l) clearly satisfy the constraints given in Eqs. ([T5|) and (fl6j) . 
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The pulse sequence shown in Figure [T] thus gives the full (not merely locally 
equivalent) CNOT solution for the 3-qubit DF subsystem. Removing the pi, 
p 2 , —pi, and 1 — p 2 gates results in a locally equivalent CNOT solution of 18 
pulses in 11 time steps for the DF subsystem, which is shorter than the locally 
equivalent CNOT pulse sequence for the DF subspace only [TO] . 

4 DFS Leakage Reduction Pulse Sequence 

Given a possibly leaked DFS qubit A and a fiducial, unleaked DFS qubit B, the 
DFS leakage reduction operator has the following specification. If A is unleaked 
(a superposition of basis states 1-4 in Table [JJ, leakage reduction leaves the 
encoded quantum number Sa.i,2 unchanged, but may possibly alter the gauge 
quantum number Sa,z- If A is leaked, leakage reduction returns A to any state 
in the unleaked subspace. The state of DFS qubit B after leakage reduction is 
unconstrained, and will generally be leaked. Since we desire an exchange-only 
leakage reduction operator, the second, fiducial DFS qubit B is required. A 
single, leaked DFS qubit has total spin S = |, which cannot be converted to an 
unleaked S — \ state by exchange operations on the single DFS qubit alone. 

Though the leakage reduction operator returns leaked states to the valid 
unleaked subspace, the final unleaked state is not necessarily a valid subsystem 
state — the reset state will not generally be a factorizable state of the encoded 
qubit and gauge qubit. After leakage reduction of a leaked state in DFS qubit 
A, the A wavefunction will be 

\ipA) = ai\ip\}\(/)gj} + axX\\i/)\)\<f>g t x) + ct Y Y x \^\)\(j>gy) + a z Z x \ip\)\cj>g,z), 

(20) 

where are complex amplitudes, \ip\) is the wavefunction of the encoded qubit, 
\4> g .^) are wavefunctions of the gauge qubit, and X\, Y\, and Z\ are single 
encoded qubit Pauli operators acting on the encoded wavefunction alone. The 
leakage reduced state thus generally appears as a superposition of the correct 
encoded state with Pauli-errors on the encoded state coupled to different gauge 
states. Standard quantum error correction will return such a superposition state 
in the unleaked subspace to a valid factorizable subsystem state. Quantum error 
correction procedures are constructed from DFS encoded gates such as the DFS 
CNOT, Hadamard, and Pauli gates. All the DFS gates have been constructed 
to act independently of the gauge state, so that quantum error correction on 
the encoded parts of Eq. ([20)1 proceeds as usual. A projective measurement in 
the course of the error correction procedure selects one of the terms in Eq. ([2"0]l . 
resulting in a valid factorizable DFS state. A similar argument holds in the case 
of fully coherent error correction. 

The leakage reduction operator in the total angular momentum basis must 
have the following form. In the spin-0 subspace the leakage reduction matrix 
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where {di j } is an arbitrary 2x2 unitary matrix. The upper left 4x4 block is the 
outer product of the identity on the encoded quantum number of DFS qubit A 
and an arbitrary unitary matrix on the encoded quantum number of DFS qubit 
B. It ensures that the encoded quantum number of DFS qubit A is unchanged, 
with no constraint on the unitary evolution of DFS qubit -B's encoded quantum 
number, in the S'tot = 0, S^.tot = gauge state. The decoupling of the fifth 
basis vector from the others ensures that no leakage from DFS qubit A occurs. 
In the spin-1 subspaces the leakage reduction matrix must be 





ei,i 


ei,a 














/i,iei,3 


/i,2ei,3 


\ 




e 2 ,i 


e2.2 














/l,l e 2,3 


/l,2C2,3 













ei,i 


ei,2 








/2,iei,3 


/2,2ei,3 













e2,i 


e2,2 








/2,ie2,3 


/2,2C2,3 







e3,i 


e3,2 














/l,l63,3 


/l,2e3,3 













e3,i 


e3,2 








/2,ie3,3 


/2,2e3,3 



















5i,i 


31,2 








31,3 
















32,1 


32,2 








32,3 


V 














33,1 


33,2 








33,3 / 



where {eij}, {fij}, and {gij} are all arbitrary unitary matrices. The ordering 
of matrix entries again corresponds to the basis vector ordering in Table OS 
Column 1 of Eq. (|2"2"|) gives the action of the leakage reduction operator on basis 
vector 6 in Tabled Basis vector 6 can be brought to basis vectors 6,7, and 10, in 
any complex linear combination: this action preserves the encoded quantum 
number of DFS qubit A and keeps A unleaked, while allowing DFS qubit B to 
change its encoded state, or even leak. Basis vector 7 (column 2) has a similar 
evolution as basis vector 6. Columns 3 and 4 must be the same as columns 
1 and 2, albeit with shifted entries; this ensures the identity is applied on the 
encoded quantum number of DFS qubit A. The constraints on columns 1-4 
and the constraints on the spin-0 matrix together ensure that the identity is 
applied on the encoded quantum number of DFS qubit A regardless of gauge 
state, assuming an unleaked DFS qubit B. 

The remaining columns of Eq. (|2"2"|) repair leaked states of DFS qubit A in 
the spin-1 subspaces. Columns 7 and 8 (basis vectors 12 and 13) bring a leaked 
DFS qubit A into basis vectors 6-11, which are all unleaked in DFS qubit A, 
leaving no components in the leaked basis vectors 12-14. Unitarity of Eq. (|22[) 
forces the entwining of the {e^} and {fi.j} matrices. Unitarity again forces 
columns 5,6, and 9 to have entries in their first six components, leaving {<?ij} 
to be an arbitrary 3x3 unitary matrix. 

Leaked states of DFS qubit A also exist in the spin-2 subspaces and must be 
reduced to unleaked states. In the spin-2 subspace the leakage reduction matrix 
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must be 



hi,i 

h 2 ,x 

h, 2 o 

k 2 ,2 

k 3 , 2 



hi, 2 \ 





\ fc 3,l 



h 2<2 

fel,3 

fc 2 , 3 

fe 3 , 3 / 



(23) 



where and {kj} are arbitrary unitary matrices. The ordering of matrix 

entries in Eq. ([23)) corresponds to the basis vector ordering in Table |U Columns 
3 and 4 ensure that basis vectors 17 and 18, which are leaked in DFS qubit A 
and unleaked in DFS qubit B, are brought to basis vectors 15 and 16, which are 
unleaked in DFS qubit A. We stress that here and in the other spin subspaces, 
DFS qubit B must be unleaked, or leakage reduction of DFS qubit A will fail. 
The spin-3 subspace of the leakage reduction operator acts only on basis vectors 
that have both DFS qubits leaked and is unconstrained aside from unitarity. 

An objective function for the leakage reduction operator is constructed from 
the constraints given in Eqs. ([2"Tj) - ([2"3l . We additionally constrain {fij} — 
{hi j}, which forces a leaked state in DFS qubit A to be reset to a factorizable 
state of encoded and gauge quantum numbers. (This constraint is not necessary, 
but enabled us to find an analytic solution.) Again, let U be the product 
of exchange gate unitaries in the total angular momentum basis, and order 
the basis vectors according to Tables HHU We define the following matrices, 
constructed from components of U : 




(24) 



(25) 



(26) 



(27) 



(28) 



(29) 



with = HlL ( ti, and 




(30) 



(31) 
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Figure 2: Thirty pulse, 20 time step, exchange gate sequence for DFS leakage reduction operator. DFS qubits A (potentially 
leaked) and B (fiducial, unleaked) are arranged as shown. Swap powers are displayed explicitly in the gates. q\ = arccos(l/3)/7T 
and 52 = arcsin(l/v / 3)/7r. 



With these definitions the leakage reduction operator objective function is 



fuvo(U) 



~(dw + dW?.(dW+dP>) 



lfft.fr _/ 2 ||, 



(32) 



where /„ is the n x n identity matrix. The objective function is zero when: 
= and each is unitary, satisfying Eq. (|2"Tj) ; H is unitary, satisfying 
Eq. ([23"]>: and = £ (2) and each is unitary, satisfying Eq. ([22} and the 

additional constraint that {fij} = {hij}. 

Figure [2] shows the best solution found by the genetic algorithm, requiring 30 
pulses in 20 time steps. Again, the solution found is analytic, with the exchange 
gate powers shown in the figure. In the total angular momentum basis the 
product of the exchange gate unitaries satisfies the forms given in Eqs. (|21[) - 
(|23|) . In Eq. (|2"Tj) the free phase satisfies e 1 ^ = 1, while the constituent unitary 
matrices, up to a common global phase, are explicitly: 



{ d i,j} 
{ e i,i> 

{/,..! 

{9i,j} 
{hij} 



\i (i + V2) (3i + V3) 

(_l)5/6 o 



V3) 



\i (i + V2) (3i + V3) 
n _2\/I 

u 3 



±V2H + V3) 

± (-t + 2^/2) (3i + V3) 35 (l - iv^) (3i + V3) 



■i + V3) 



1 






2 






7 




1 


972 


243 


18 


140 


80i\/2 


2\/5 


243 


243 


9 



-| (-Z + 72) (-t + V3) 
\ 



t\/2 
3 



4\/5 
9 

_ _L 

9 



( 



J_ , n/2 
12 3 



_7 |y! 

12 3 





\ 


-1 / 



(33) 
(34) 
,(35) 

(36) 

(37) 
(38) 



5 Genetic Algorithm 

The DFS CNOT and leakage reduction operator gate sequences were con- 
structed using the genetic method to minimze objective functions associated 
with the desired CNOT and leakage reduction operator block matrix forms. 
Although simple genetic algorithms can be quite slow, finding sequences of the 
lengths needed for the CNOT and leakage reduction operators (fewer than 40 
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gates) can be accomplished using a rather naive programming approach in a 
few weeks of evolution on an Apple XServe computer. (We also attempted opti- 
mization by simulated annealing, but the gate sequences we present here for the 
CNOT and leakage reduction operator were both found by genetic methods.) 
It helps that the time consuming part of the computation is embarrassingly 
parallel. 

Our genetic approach iterates random changes on some population of gate 
configurations, each given by a list of powers of swap between adjacent phys- 
ical qubits. At each iteration ("generation") we augment the population with 
"mutations" and "mating" , and then apply "natural selection" that favors those 
configurations with low value of the objective function. Wc augment the ob- 
jective functions given in Eqs. (fT7| and (|32|) by adding a "gate penalty" term 
which is simply the length of the sequence times some positive constant, the 
gate penalty parameter. The gate penalty parameter is adjusted as the evo- 
lution progresses to give a reasonable tradeoff between lowering the objective 
function and lowering the gate count. 

The code was written in Mathematica Version 7.0 [18] . using the built-in 
parallel processing capability. The two key Mathematica functions used were 
FindMinimum and RandomSample. FindMinimum searches for a local minimum 
of a function using a variety of local descent methods; RandomSample gives a 
pseudorandom sampling of a list with optionally weighted probability. Further 
description of these functions may be found in the Mathematica documentation 

5.1 Mutations 

We have programmed five different types of mutations, chosen at random. We 
have varied the probabilities of each mutation, but for our CNOT and leakage 
reduction operator solutions, the probabilities are as shown: 

Ref ineOne (8°/ ) : randomly chooses one gate of a configuration and uses 
FindMinimum to optimize its power of swap, starting with the current 
value. 

Ref ineTwo (87.) : randomly chooses an adjacent pair of gates and simultane- 
ously minimizes the objective with respect to both powers of swap. 

Ref ineAll (4°/ ) : runs FindMinimum on all powers in the configuration. The 
time required to do this grows rapidly with the number of variables and is 
much larger than that of the previous two mutations. The Mathematica 
AccuracyGoal and PrecisionGoal Options of FindMinimum are set to 
modest values (< 5) so that this operation runs in a tolerable length of 
time. 

InsertGate (40%) : inserts a gate at a random location with a random power, 
followed by Ref ineOne. 
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DeleteGate (407 ) : deletes a gate at a random location, followed by Ref ineTwo 
on the gates that were adjacent to the one removed. If the removed gate 
was the first or last, the Ref ineOne is used. 

5.2 Mating 

The initial population (unmutated) is divided into pairs. We take the begin- 
ning of one gate configuration and append the end of the other, and then run 
Ref ineAll. The length of the resultant "offspring" is chosen to be the mini- 
mum of that of the parent configurations and the number of gates from each 
parent is chosen at random. Because of Ref ineAll, this is a time consuming 
operation - but without the Ref ineAll very few of the offspring would survive 
natural selection. The genetic algorithm can be used without mating, but the 
diversity of the population suffers. 

5.3 Natural Selection 

Surviving gate configurations are selected using the Mathematica function 
RandomSample. The sampling weight, which is proportional to the probabil- 
ity of a member of the population appearing in the sample, is given by the 
inverse of the objective function augmented by the gate penalty term. Config- 
urations with a small value of the objective function and low gate count are 
most likely to survive. The population size is specified by a parameter under 
user control. After the survivors are determined, gates between the same pair 
of qubits that are not separated by any noncommuting gates are combined into 
one gate, by simply adding their powers. 

6 Conclusion 

Using a genetic algorithm, we have found new, exact exchange gate pulse se- 
quences for the encoded CNOT operation on two 3-qubit DFS qubits, as well 
as for the DFS leakage reduction operator. Our CNOT solution gives the most 
efficient sequence yet found for either the 3-qubit subspace or subsystem, while 
our leakage reduction solution is the first explicit sequence for any exchange- 
only subspace or subsystem encoding. We have found a wide range of other 
solutions for the DFS CNOT and leakage reduction operator, but the sequences 
presented here are the most efficient thus far. Other solutions that place ad- 
ditional constraints on the evolution of the leaked states may have superior 
error propagation properties when combined with fault tolerant error correction 
procedures; further investigation of such sequences is underway. 
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